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Delayed dynamic triggering of earthquakes: Evidences from a sta- 
tistical model of seismicity 
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Abstract. - I study a recently proposed statistical model of earthquake dynamics that incorpo- 
rates aging as a fundamental ingredient. The model is known to generate earthquake sequences 
that quantitatively reproduce the spatial and temporal clustering of events observed in actual 
seismic patterns. The aim of the present work is to investigate if this model can give support to 
the empirical evidence that earthquakes can be triggered by transient small perturbations, partic- 
ularly by the passing of seismic waves originated in events occurring in far geographical locations. 
The effect of seismic waves is incorporated into the model by assuming that they produce instan- 
taneous small modifications in the dynamical state of the system at the time they are applied. 
This change in the dynamical state has two main effects. On one side, it induces earthquakes that 
occur right at the application of the perturbation. These are called immediate events. On the 
other side, after the application of the perturbation there is a delayed effect: the seismic activity 
increases abruptly after the perturbation, then falls down below the level of background activity, 
and eventually recovers to the background value. The time scale of these variations depends on 
the internal dynamics of the system, and is totally independent of the duration of the perturba- 
tion. The number of delayed events in excess of the background activity is typically observed to be 
around a factor of twenty larger than the number of immediate events. The origin of the enhanced 
activity period following the perturbation is associated to the existence of aging relaxation, and 
it does not occur if relaxation is absent. These findings give support to the experimental evidence 
that earthquake can be remotely triggered by small transient perturbations as those produced by 
seismic waves. 



Introduction. — Some earthquakes are triggered by 
previous earthquakes. The best known example is the oc- 
currence of aftershocks, which are triggered by a previous 
large event, called the main shock. Aftershocks originate 
in the stress redistribution that the main shock produces 
around the rupture region Q] . This is why aftershocks oc- 
cur typically inside or near the rupture region of the main 
shock 2 . Aftershocks are an example of 'static trigger- 
ing', since the stress redistribution that generates them is 
permanent (at least until stresses are modified by future 
events) . 

On the other hand, 'dynamic triggering' refers to the 
possibility that some earthquakes are triggered by tran- 
sient perturbations of the stress field, most remarkably 
by the passing of seismic waves produces by other earth- 
quake [3J. Observational evidence of dynamic triggering 



has been accumulated since the nineties. For instance, 
the 1992 Landers earthquake triggered seismic activity as 
far as Yellowstone National Park [4] , and the 2002 Denali 
earthquake triggered seismicity in southeastern California 
[5J. Following these early observations, many evidences of 
dynamic triggering at large spatial distances have been ac- 
cumulated (see [3] for a detailed set of references), in such 
a way that by now dynamic triggering of earthquakes is 
widely accepted to occur. Recently, evidences of dynamic 
triggering have been found also at the laboratory scale [5] . 
Dynamically triggered events may occur right at the pas- 
sage of the seismic waves [4,5,7. . This is a case of 'instan- 
taneous' triggering. However, one of the most intriguing 
characteristics of dynamic triggering is that in most cases 
the triggered events occur much later (hours or days) than 
the perturbation. This kind of 'delayed' dynamic trigger- 



p-1 



E. A. Jagla 



ing is revealed by a comparison of the average seismicity 
before and after the passage of the seismic waves, and is 
undoubtedly documented for instance following the 1992 
Landers event, which was seen to activate seismicity in a 
large region of western North America [1] . 

The qualitatively most accepted explanation of dynami- 
cal triggering involves two steps. The first one corresponds 
to assume that the passage of seismic waves produces some 
kind of damage, or modification in the state of the fault, 
that contributes to 'advancing the clock' of its seismic cy- 
cle [8]. The kind of physical modification that seismic 
waves can produce may involve for instance the increase 
in pore fluid pressure, which leads to an effective decrease 
in effective normal stresses, favoring a weakening of the 
fault [3]. Another possibility for the physical effect of the 
passing seismic wave is the acceleration of subcritical crack 
growth 0: a sudden, although transient, transfer of stress 
to the rock will cause the size of the cracks to increase some 
amount, moving the state of the fault towards a situation 
of instability. The second step in the dynamical triggering 
process is the completeness of the seismic cycle of the fault, 
leading to failure at an earlier time that it would have been 
in the absence of the perturbation. On a theoretical per- 
spective, the process of advancing the clock of the fault 
and its eventual failure have been described in terms of 
the phenomenological rate-and-state equations modeling 
fault friction [HUD], and the results obtained support the 
previous schematic description. 

Beyond this rather qualitative understanding of delayed 
dynamic triggering, it would be important to have a model 
in which the described scenario can be verified in its de- 
tails, without introducing special hypothesis related to dy- 
namic triggering itself. Statistical models, for instance 
spring-block, or similar kind of models, have in the past 
provided a solid basis for the understanding of basic fea- 
tures of seismic phenomena [TTtfT4] . For instance, it was 
found that the broad distribution of earthquakes sizes ob- 
served, and accounted for by the Gutenberg-Richter law 
[2], can appear as an emergent collective behavior in a 
system consisting of a large number of elements having 
all the same typical size. Recently [THUS], ^ has been 
shown how the inclusion of internal relaxation, or aging 
in these models, expands greatly their possibilities, giving 
rise to sequences of events that display typical spatial and 
temporal correlations compatible with actual sequences, 
and also showing that relaxation can be considered as a 
physical ingredient that generates friction properties com- 
parable with those experimentally observed. 

The aim of the present work is to study the possibility 
of dynamic triggering of earthquakes in this kind of sta- 
tistical model with relaxation. As in the interpretation of 
field results, there is here also the problem to define on 
safe grounds what the actual effect of the perturbation is 
on the state of the system. This will necessarily imply 
some ad hoc hypothesis to be added to the definition of 
the model. I will assume some kind of random damage 
(explained in detail below) caused by the passage of the 



seismic wave, and take this damaged configuration to be 
the starting point of a simulation that is compared with a 
situation in which the perturbation is not included. The 
main results that I obtain are the following: The pertur- 
bation triggers some immediate events (occurring right at 
the time when the perturbation is applied). In addition, 
during the temporal evolution after the application of the 
perturbation, a variation of the average seismicity rate in 
the system is clearly detected as compared to the case 
without perturbation. Seismic activity is at a maximum 
right after the perturbation, decreases in time becoming 
lower than the background reference value, and then re- 
covers to the reference value. The typical time scale of this 
process is comparable to the time scale of the aftershock 
cascades of the largest events observed in the system with- 
out perturbation, and this time scale has nothing to do 
with the time during which perturbation is applied. The 
time-integrated temporal density of events does not differ 
significantly between the cases with and without pertur- 
bation, and in this sense, the perturbation can be though 
of as advancing in time the occurrence of some events in 
the system, instead of generating new ones. By counting 
the events in the peak of over abundance with respect to 
background after the perturbation, this number is seen to 
be around a factor twenty larger than the number of im- 
mediate events triggered by the perturbation, pointing to 
the fact that delayed events will be much more significant 
than immediate events during dynamic triggering. The 
perturbation does not alter the distribution of events in 
magnitude from the case without perturbation. A small 
perturbation produces a rather large effect in the system 
in the form of delayed events: a perturbation that affects 
1% of the sites of the system produces eventually delayed 
events involving about 20% of all sites. 

In addition, in the range of parameters studied, the vari- 
ations of the seismicity observed around the background 
value are proportional to the intensity of the perturbation, 
namely the effect of the perturbation can be reasonably 
described as the superposition of the effects produced by 
the damage produced on individual sites by the perturba- 
tion. Before presenting the results of the simulations, for 
completeness in the next section the basic assumptions of 
the model are explained. 

Brief recalling of the model, and the damage as- 
sumption for dynamical triggering. I give here the 
basic definitions of the model, as it was presented in |16j . 
An alternative justification from a more physical starting 
point with essentially the same results can be seen in [13] . 
The model can be considered to be an elaboration on the 
Olami Feder and Christensen (OFC) model to describe a 
single planar fault [13] . In the OFC model a real variable 
fi is defined on every node of a two dimensional square 
lattice, represented here by the index i. The value of fi is 
interpreted as the local value of the friction force between 
a sliding block and a rigid underlying surface. Tectonic 
loading at some velocity V is applied by increasing the 
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values of all fa linearly in time, i.e. dfi/dt = V. Each 
time one of the fa reaches a threshold value of 1, a local 
instability occurs and an unload process occurs. The cor- 
responding fa is set to zero, and the neighbor sites j are 
loaded according to fj — > fj +afa, where a is a parameter 
between and 1/4. This can produce some neighbor sites 
to overpass the threshold, and so a cascade can occur, un- 
til all fa < 1. This cascade is an earthquake, or an event 
in the system. 

The modifications introduced in [TB] to the OFC model 
to obtain realistic features of earthquakes are the follow- 
ing: 1- thresholds are not uniform but they are site de- 
pendent, noted by and chosen from a Gaussian distri- 
bution with unitary mean and deviation a (in this paper 
I will use a = 0.3). Each time a site reaches the thresh- 
old and is downloaded, /, is decreased by 1, an amount a 
is uploaded onto each neighbor and the local threshold is 
chosen anew. The magnitude of each event in the system 
is calculated as M — 2/31og 10 5, where the 'seismic mo- 
ment' S is the total number of topples that composed the 
complete avalanche. 2- The time evolution of fa between 
events is not simply governed by tectonic loading, there is 
now a relaxational term that tends to make the values of fa 
spatially uniform. The equation that describes temporal 
evolution of fa is 

y±=V + R(V 2 V 2 f) i (1) 
dt 

where V 2 is the discrete Laplacian on the square lattice. 
This is one of the possible relaxation terms discussed in 
[IT)] . A physical motivation of its form can be found in 
[T5] . Note that the alternative of using a single Laplacian 
instead of the double one in eq. ([TJ produces only minor 
differences, as explained in |16j . The relaxation parameter 
R controls the intensity of the relaxation process. In ad- 
dition, periodic boundary conditions are used. The model 
just described generates temporal sequences of events that 
display realistic aftershocks, produces a GR law with a 
correct decaying exponent without tuning of parameters, 
and displays realistic features of the friction process, as 
for instance the phenomenon of velocity weakening. 

On top of this model, I want to study the effect of the 
perturbation produced by the passage of seismic waves. 
These seismic waves are supposed to be generated at some 
remote location that is not included in the modeling. In 
particular, they do not originate in the earthquakes that 
are observed to occur in the model. An additional hypoth- 
esis is necessary here, since we do not really know what the 
effect of the passage of seismic waves (or any other pertur- 
bation) is. But is not totally unreasonable to assume that 
the effect, whatever the actual physical mechanism is, can 
be described as a (small) change in the dynamical state of 
the system. I consider two possible modifications gener- 
ated by the perturbation. They will be called model A and 
model B. In model A, at the time of the perturbation, a 
fraction e of the sites of the system are chosen at random, 
and they are unloaded by the normal rules in the model. 



This may be interpreted as if there is a uniform probabil- 
ity e of failure caused by the passage of the seismic wave. 
In the second realization (model B), the sites that fail are 
those for which /*" — fa is lower than some threshold value 
A/. This may be interpreted as the seismic wave produc- 
ing the oscillation of the values of fa with some amplitude 
A/. In this case, the intensity of the perturbation can be 
characterized also by the fraction e of sites that became 
destabilized because of the perturbation. 

Results. — Before going into the detailed results, it 
is convenient to emphasize that the delayed dynamic trig- 
gering to be observed in the model is associated to the 
existence of relaxation (i.e., a non zero value of R in eq. 
(JTJ). In fact, in the case in which R = 0, the pertur- 
bation causes no variation of the average activity, except 
for the existence of some immediately triggered events (in 
an amount similar to those observed with relaxation), in 
particular, delayed triggered events are not detected for 
R=0. 

Dynamic triggering can be characterized as the differ- 
ence in the seismic activity of the system when a per- 
turbation is applied, compared to the case in which no 
perturbation was applied. In fig. [TJa) we see an example 
of the temporal activity in a system of 256 x 256 sites in 
which a perturbation was applied at t = with an inten- 
sity e = 0.01 (i.e., one percent of the sites were instanta- 
neously destabilized by the perturbation) . It is clear from 
this figure that due to the intrinsic fluctuations in seis- 
mic activity, we cannot extract any reliable conclusion on 
the influence of the perturbation on the dynamics. A pos- 
sibility to circumvent this problem would be to increase 
the system size, so the fluctuations in the seismic activ- 
ity average out. But this would imply to work with pro- 
hibitively large (in terms of computational time) system 
sizes. A more reasonable possibility is to consider many 
realizations in which perturbation is applied, and seismic 
response is presented as the average value over realiza- 
tions. In fig. [T](b) and (c) we see the result of accumu- 
lating many realizations for the two kind of perturbations 
discussed, namely model A and model B. Now a clear dif- 
ference between activity with and without perturbation 
emerges, which is the phenomenon of dynamic triggering. 
To quantify in more detail distributions of events like those 
in fig. [Ub-c), I will present results mainly for two quan- 
tities: temporal density of events N(t), i.e., the number 
of events occurring per unit of time, and temporal density 
of stress drop er(t), i.e., the seismic moment of all events 
occurring per unit of time. For both quantities, a lower 
size cutoff So will be indicated, i.e., only events with seis- 
mic moment S > Sq will be considered. By definition, 
both N(t) and a(t) are constant in time in the case of no 
perturbation. After a perturbation is applied, there will 
be changes in these quantities, which recover their asymp- 
totic values at large times. 

In fig. [2]we see the results for N(t) and a(t) for accumu- 
lated sequences like those in fig. [TJb-c). The results are 
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Fig. 1: (a) A single realization of the events in a system of size 
256 x 256 (R/V = 50, a = 0.24). At t = a perturbation 
(of type A) was applied, which destabilized a fraction e — 0.01 
of the sites in the system. This produced at the same t = 
cascades that involved a fraction of about 0.03 of all the sites 
in the system. These immediate events are not shown in the 
figure, (b) The result of accumulating 50 realizations like the 
one in panel (a), (c) The equivalent results accumulated during 
50 realizations, now for a type B perturbation, with the same 
value of e = 0.01. In (b) and (c) the overabundance of events 
following the perturbation is clearly discernible, as well as a 
broader region of depleted activity at longer times. 



normalized to the background activity, and curves are pre- 
sented for different values of the low size cut off So used in 
counting events. The perturbation applied at t — brings 
about two main effects. First of all, there is a number 
of events that occur exactly at the time of the perturba- 
tion. They can be understood as the immediate cascades 
triggered by the sites that become unstable by the appli- 
cation of the perturbation. In addition, once the pertur- 
bation has disappeared and as time increases, the seismic 
activity is first enhanced, then it decreases, becoming even 
lower than the unperturbed value, to finally rejoin the un- 
perturbed value at very long times. The events that form 
the peak above the unperturbed value will be referred to 
as the delayed dynamically triggered events. Their exis- 
tence is in fact, the single most important result that I 
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Fig. 2: Seismic activity as a function of time, averaged over 
about 400 realizations, following a perturbation occurring at 
t = 0. Results for different values of the lower size cutoff So 
are indicated, to show the independence on this cutoff. The 
results are normalized to the corresponding background activ- 
ity. Immediate events are not displayed here, but they amount 
systematically to about 5% of the number of delayed events 
(i.e., the peak above the background activity values) 



present in this paper. Within the precision of the curves 
in fig. [5J the results are independent of the low size cutoff 
So used to count the events, both for the immediate events 
and for the delayed ones. This is a general result that was 
observed in all simulations: dynamic triggering produces 
temporal variations of the seismic activity, but the distri- 
bution of events in magnitude continuous to be an intrinsic 
property of the model, independent of the perturbation. 
The time integrated values of both N t and a(t) are (again 
within the available precision) equal to the values without 
perturbation. For a(t) this is not a surprise, since on very 
large time intervals, the total stress drop can be calculated 
as the applied external stress on the system, and this is in- 
dependent on the existence or not of a perturbation of the 
kind we are considering. For the case of N(t) the fact was 
not obvious from the beginning. The fact that the number 
of events in the enhanced and depleted regions in fig. [2] 
are almost equal, comes out as a non-trivial results, which 
is actually related to the fact already mentioned that the 
distribution in magnitude of events is not altered. 

The fact that the perturbation does not change the to- 
tal number of events allows a possible interpretation of the 
dynamical triggering in the following form: It can be said 
that dynamical triggering proceeds by making some events 
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that were 'scheduled' to occur later in time, to appear ear- 
lier, because of the perturbation. While this interpreta- 
tion is possible in principle, I think it is a bit dangerous 
to take it too seriously, as there is no clear comparison 
between the activity observed in a single realization of the 
dynamics with and without the perturbation. 

In comparing the results for model A and B in fig. [51 
we see that the results are qualitatively similar in both 
cases. A quantitative difference is that the time scale of 
the process is longer for model A, and shorter for model 
B. The reason of this behavior can be attributed to the 
fact that model A destabilizes a fraction of sites at ran- 
dom, whereas in model B the sites that are already close to 
instability are those effectively destabilized by the pertur- 
bation. This produces that for model B a lower amount of 
tectonic loading {i.e. a shorter time) is necessary to erase 
the effect of the perturbation as compared to model A. 

In addition, in a broad range of perturbation intensity 
studied, the effect of the perturbation is simply propor- 
tional to the perturbation intensity. In fact, results ob- 
tained for different values of the intensity e of the per- 
turbation (for e between 0.001 and 0.02) show that the 
amount by which the results differ from those in the ab- 
sence of perturbation scales linearly with the perturbation 
intensity. In particular, this implies that the time scale in 
which the delayed events are observed is independent of 
the perturbation intensity. This linearity of the effect with 
the intensity of the perturbation, allows to make a concise 
statement on the number of events triggered by a pertur- 
bation. In fact, the number of immediate events observed, 
amounts typically to roughly 5% of the number of dynam- 
ically triggered delayed events. This result is independent 
of the value of So used as a lower size cutoff. Another way 
of stating this is to say that the eventual effect that a per- 
turbation has onto the system is about a factor of 20 larger 
than what would be naively expected, i.e., than the effect 
observed immediately after the perturbation disappeared. 
This is an important result, because it points to the par- 
ticular sensibility of the system to the perturbation, and 
indicates that (as it actually occurs in real seismicity) the 
largest amount of triggered events occur after a certain 
delay, and only a small fraction is instantaneous with the 
perturbation. 

Additional insight into the process of dynamical trigger- 
ing is obtained by studying the response of systems with 
different intensity of relaxation R. First of all, it is clear 
from eq. (HJ that a simultaneous change in V and R by 
the same factor is equivalent to a change in the time scale 
by the same factor, i.e., upon such a change in V and 
R, results as those in fig. [5] would look exactly similar, 
with only the time axis rescaled. However, an indepen- 
dent change of R and V produces characteristic effects in 
the response of the system. For instance, in fig. [3] we see 
curves for different values of R, keeping V as fixed. It 
can be seen that the temporal duration of the peak of dy- 
namically triggered events is roughly proportional to 1/R, 
but its area (i.e., the total number of dynamically trig- 
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Fig. 3: Results for different values of the relaxation parameter 
R in eq. |T]) (V = 1 and e = 0.01 and So = 1 in all cases, results 
for model A). In (b) the axis are scaled according to the value 
of R, and the data are presented in logarithmic scale. It can 
be seen that the peak of delayed dynamically triggered events 
scales with R in such a way that its area is rather independent 
of R, particularly when R becomes large, and the time span 
of these events is governed by 1/R. Instead, the time required 
for the system to recover the background activity value is in all 
cases about 0.1, independently of R. In (b), a temporal decay 
with an exponent p = 1.5 is included, for comparison. 



gered events) remains the same. On the other hand, the 
time that takes the system to recover the background ac- 
tivity is independent of R (actually, it is proportional to 
X/V). The qualitative explanation of this behavior is the 
following: dynamically triggered events are produced by 
the relaxation mechanism (the R term in eq. (HJ), and 
they would occur even in the case in which V = 0. In 
this limit, it is clear that the effective time scale is in- 
versely proportional to R, and this justifies the scaling of 
the peak of dynamically triggered events with R. On the 
other hand, once all dynamically triggered events have oc- 
curred, the system is in a state in which there has been 
an excess of seismic moment released. This deficit has 
to be eventually compensated by the loading mechanism, 
which is accounted for by the V term. Thus we see that 
the temporal extent of regions of enhanced and depleted 
seismic activity are controlled by the independent param- 
eters R and V, respectively. In addition, if the ratio R/V 
becomes sufficiently large, the enhancement of seismic ac- 
tivity right after the perturbation will be much larger than 
the largest activity depletion, which becomes a very broad 
in time and quite shallow. This limit is compatible with 
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actual seismic observations, as in this case there is no clear 
evidence of a period of depleted activity following the pe- 
riod of triggered activity. 

There is some experimental evidence that the form of 
the time decay of excess activity caused by dynamic trig- 
gering follows approximately the same Omori law that fol- 
low the usual aftershocks caused by static triggering. In 
the present model this is also the case. The decay of N(t) 
may correspond to an Omori decay (fig. Eth)). However, 
for the ratios between R and V investigated, this decaying 
is masked at large times by the region of depleted activity, 
so at this point it can be only stated that the time decay 
observed of the peak of dynamically triggered events is not 
incompatible with the Omori law. A more definite state- 
ment could be obtained by analyzing cases with larger 
values of the ratio R/V. 

Conclusions. — In the present paper I have shown 
that the same mechanism that was used in Refs. [15|I16] to 
obtain realistic sequences of earthquakes, is able to provide 
a consistent description for the phenomenon of dynamic 
earthquake triggering. A modification of the dynamical 
state of the model at some time t = 0, assumed to be 
originated in the passage of seismic waves of some remote 
event, it was shown to produce in the model the follow- 
ing set of observations: A number of "immediate" events 
at the time of the perturbation are triggered. An addi- 
tional number of events (in excess of the normal activity) 
accumulate in the time after the perturbation has already 
vanished. These delayed dynamically triggered events are 
typically a factor of 20 larger than the number of immedi- 
ate events. The time scale of delayed events is controlled 
by the relaxation parameter R in the model, and in par- 
ticular is independent of the duration of the perturbation 
(here it was assumed that perturbation occur during a 
vanishingly small time interval) . The perturbation has no 
effect on the magnitude distribution of events. The effect 
is proportional to the intensity of the perturbation, in a 
broad range of this intensity. 

These findings give support to the somewhat debated 
phenomenon of dynamic triggering. The central point 
that the model captures, is the way in which a transient 
and small perturbation is able to trigger events that occur 
much later in time, once the perturbation has disappeared. 
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